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00 '. Abstract 

7— I ■ We demonstrate that the relative size of the CP-even and CP-odd cou- 

plings of a Higgs boson to ti can be determined at the LHC at a level that is 
very useful in discriminating between models, provided that a Higgs signal 
C*~) ' can be observed at a reasonable S/^fB level in the tth (with h — > 77 or 

bb) channels. In particular, the CP-even nature of a SM-like Higgs boson's 
\ coupling to ti can be confirmed with very substantial statistical significance 

Q\ I using the tth (h — ► 77) production/decay mode. The key to achieving good 

■ discrimination between models is to make full use of the difference in the 

Oh, final state distributions for different ratios of the CP-even to CP-odd cou- 

Q_i! plings. For our analysis, we employ the very convenient and simple optimal 

analysis procedures derived previously. 

4^ 



X : 1 Introduction 



It will be very important to directly determine the CP nature of any Higgs boson 
(h) that is discovered. Since it is almost certain that a substantial number of Higgs 
boson events will first be available at the LHC, it will be very desirable to use LHC 
data to verify the CP nature of any observed Higgs without waiting for an e + e~ 
collider. The only procedure with a reasonable level of viability proposed to date 
is to employ weighted integrals of the pp — > tthX final state JTJ. These were shown 
to provide significant ability to establish that a SM-like Higgs boson has a purely 
CP-even ti coupling. Subsequently, procedures for optimizing the extraction of the 
relative magnitude of different components of a cross section were developed and 
their power demonstrated when applied to determining the relative magnitude of 
the CP-even and CP-odd tth couplings in the process e + e _ — > tth 0. Here, we 
apply these same techniques to the tihX final state at the LHC and demonstrate 
substantial improvement as compared to the weighted integral technique of Ref . [ffl] . 
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In general, a very meaningful determination of the relative magnitude of the CP- 
even and CP-odd tt couplings of the h will be possible provided only that a signal 
with reasonable statistical significance is present. The key to the optimal analysis 
techniques is to take full advantage of all the information available in the cross 
section as a function of the kinematical variables. 

For our analysis, it is important that a relatively pure sample of ttFX events 
can be isolated for one or more final states F into which the Higgs decays. Here, 
we will consider the F = 77 and F = bb final states originally discussed in Refs. || 
and UJ, respectively. The events will consist of tihX (h — > F) signal events and 
background events. (For a light h, our focus here, interference between signal and 
background amplitudes can be neglected.) The pp — > ttFX cross section takes the 
form 

E(0) = g = C [a 2 fM + b 2 f b {4>)] + M0) , (1) 

where denotes the final state phase space configuration, C is a coefficient deter- 
mined by experimental efficiencies, BF(h —* F) etc., and a and b are the CP-even 
and CP-odd Higgs couplings appearing in the Feynman rule, 

tth : -t(a + iby 5 )t-^- . (2) 
2m w 

(The SM Higgs boson has a = 1 and b = 0. A purely CP-odd Higgs boson has 
a = and b 7^ 0.) We note that there is no ab cross term in Eq. (P since we do not 
consider any observable sensitive to the spin directions of the top quarks. Further, 
since both the a 2 and b 2 terms are of a CP-conserving nature, the functions f a and 
fi, do not change when the t and t momenta are interchanged. 

In Eq. (jl]), f a ,b(4>) are, in principle, precisely known functions of and have well- 
determined relative normalization. The only possible uncertainties in the f a ^ arise 
through their dependence on quark/gluon distribution functions; we will assume 
that the relevant distribution functions will have been determined with sufficient 
precision, using other LHC data and data from other accelerators, that these un- 
certainties can be neglected. We also assume that the background contribution 
/b(0) will be precisely determined using experimental data outside but near the 
Higgs mass peak in the final state F. This precise measurement will include an 
accurate determination of both the functional form and the normalization, includ- 
ing all experimental efficiencies. Statistical fluctuations in the background within 
the accepted Higgs mass bin in the state F are automatically incorporated in the 
procedure we will employ. 

2 The Optimal Analysis Procedure and Results 

As in Ref. 0, we wish to discriminate between models purely on the basis of 
the difference between the phase space distribution dependence of f a and ft, not 
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relying on knowledge of C. To this end, it will be convenient to define the rescaled 
functions f a ,b,B{4>) = fa,b,B{4>) I 'I a ,b,B where I aAB = // a ,6,B(0)# so that f f a ,b,B = 
1. (The 4> integral will be restricted by appropriate cuts.) Defining r = I\,/I a and 
a = a 2 /(a 2 + rb 2 ), we may rewrite Eq. (U) as 

S(0) = 0-5 [a/„(0) + / b (0)] + o-b/b(0) , (3) 

where f a (<f>) = [/ a (0) - / 6 (0)] (J f a (<f>)d<f> = 0) and a 5 = C(a 2 / a + b 2 I b ) is the 
effective integrated signal cross section. The background cross section is ob = Ib- 
Clearly, the ratio a is independent of the overall signal normalization (although 
errors in a will certainly depend on the signal and background rates). We now use 
the optimal analysis technique of Ref. [@] to provide a determination of a with the 
smallest possible error (in the Gaussian statistics approximation). 

In general, for £(</>) = J2iCifi(4>) + g{4>) (where g{4>) and the fi(<j>) are known 
functions, including normalization) the determination of the unknown signal coef- 
ficients q with smallest possible statistical error is given by @ 

c i = y £ , where M lk = f lMM^. d(f) and j fc = f f k ^) d(f) . ( 4 ) 

The covariance matrix for the q is 

M~ l a T 

= (AcAcj) = = M^/L eS , (5) 

where <jt = J ^d<P is the integrated cross section and N = L eS a T is the total 
number of events, with L e g being the luminosity times efficiency This result is 
the optimal one regardless of the relative magnitudes of the different contributions 
to In the Gaussian statistics limit, it is equivalent to determining the Cj 

by maximizing the likelihood of the fit to the full <ft distribution of all the events. 
The increase in errors due to statistical fluctuations in the presence of background 
are implicit in the background contribution to £(</>) appearing in Eq. (|J), which 
implies larger M^ 1 entries in Eq. (|5|). From the result of Eq. (|5|), the x 2 in the Q 
parameter space is then computed as 

X 2 = J2( c 'i - c i) V ij l ^j ~ c j) = Y,( c 'i - CijLegMijic'j - Cj) , (6) 

where, for our theoretical analyses, the c, are the input model values for which 14,- 
is computed. [| 

The above procedures are not altered if cuts are imposed on the kinematical 
phase space over which one integrates; one simply restricts all <fi integrals to the 



^^Once experimental data is available, one would compute My- using appropriate binning and 
actual data for £(</>). 
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accepted region. If a subset, 0, of the kinematical variables cannot be determined, 
then the optimal technique can be applied using the variables, 0, that can be 
observed and the functions /j(0) = / fi{(p)d(j). 

We apply the above general technique to the specific case in hand by identifying 
the signal function fi = asf a and the known function g(<f>) = asfbify + 0b/b(^)- 
(The normalizations are fixed by the experimental measurements of the signal and 
background rates, which include all efficiencies etc.) The \ 2 associated with a 
model yielding a value a' that differs from the value a of the input model is then 
given by the following simple result: 

,2 /„/ „\2 & f ii fai'P) _ I I „\2^ 2 \Tjt„ C/DM2 



3 J ! [<*/«(#) + fb{<f>)] + /b(0) 3 

(7) 

where the expected signal rate S is that estimated (and eventually measured) for 
the input (true) model. It is important to note that H(a, S/B) depends implicitly 
on the channel F being considered. Defining the discrimination power, D = ^/x 1 , 
we have 

D = \ot -a\-;=H(a,S/B) (8) 
v B 



Note that in the limit of small S/B, 
H(a, S/B) -> H 



In applying the above to the pp — > tiFX process, we only consider the kine- 
matical variables associated with the tiF portion of the final state, including its 
(where s is the invariant mass-squared for the ttF combination) and overall 
rapidity. Any variables associated with the internal phase space of X are inclu- 
sively summed over as usual (and fall into the category noted above). Also, 
distributions for the decay products of the t-quark and the Higgs boson are not 
considered. 

In order to establish that we are, in fact, dealing with the ttF final state and 
completely determine the kinematical point <fi in phase space, we must be able 
to reconstruct and identify the t and t. However, identification of t vs. t is not 
actually necessary since, as noted earlier, our functions f a and /& do not change 
if the t and t momenta are interchanged. Thus, in principle we can use either the 
final state mode in which both the t and t decay completely hadronically or that 
in which one has a leptonic W decay and the other a hadronic W decay.[| In this 
paper, we choose to focus on the latter (mixed hadronic-leptonic) mode as being 
possibly somewhat cleaner. However, we anticipate that the purely hadronic final 
state could be used to further increase our statistics. 



2 The final state mode in which the Ws, from both the t and the t decay leptonically does not 
allow full determination of the kinematics of the final state. 
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In the final state mode in which one t decays leptonically and the other hadron- 
ically, there is a potential two-fold ambiguity when determining the longitudinal 
momentum of the missing neutrino by requiring that the leptonically-decaying W 
be on-shell. However, the wrong solution gives a reconstructed top-quark mass 
that differs from the known value, m t . If we assume that energy resolutions and 
such will allow reconstruction of m t to within ±10 GeV for the correct solution, 
then we can identify the correct solution in those events for which the incorrect 
solution gives a reconstructed top-quark mass that differs by more than 10 GeV 
from m t . About 70% of the events that obey the global cuts specified below satisfy 
this criterion. We accept only such events. The global cuts that we impose on the 
tth final state are such that there should be reasonable efficiency for the required 
reconstructions. We require \y t th\ < 4. A cut of pt > m h /4: is applied to the 7's 
or 6's in the F final state in order to reduce the size of the background. These are 
the same cuts as employed in Ref. [I]. 

We first give results for a Higgs mass of m h = 100 GeV, so as to allow direct 
comparison with Ref. M. (We will discuss later the sensitivity of our results to the 
\y\ cut and to m^.) For the above cuts and this mass, we find r = 0.288. (Note 
that r is independent of the final state channel F.) 

As already noted, we will consider the F = 77 and F = bb decay modes for 
the h. The F = 77 mode will only be considered in the case where the input 
model corresponds to a SM-like h, since it is only in this case that BF(h — > 77) 
is likely to be substantial. Following Ref. @, the LHC ATLAS § and CMS @ 
collaborations performed detailed simulations of this channel as part of their tech- 
nical design reports. Based on their results, and assuming an eventual integrated 
luminosity of L = 600 fb~ (summing over the two detectors), Ref. M concluded 
that one could expect S ~ 130 and B ~ 21 in the F = 77 mode for a SM-like h 
with mass of order 100 GeV. For F = bb, the best global cuts and associated levels 
of signal and background have not yet been firmly established. We will adopt the 
procedure of Ref. of examining results for several different extremes of S/B, 
assuming a certain level of statistical significance S/ \f~B. In this way, we avoid 
making any specific assumption about the branching ratio for h — > bb decay, which 
is, in any case, absorbed into the overall normalization factor C of Eq. (|IJ). For 
both F = 77 and F = bb, we will follow Ref. [0] and model the background cross 
section shape /#(</>) [see Eq. ([[])] using the irreducible pp — ► it'yyX and ttbbX 
background processes, respectively. In the case of F — 77, the ATLAS and CMS 
studies confirm the theoretical claims || that the irreducible background should, 
indeed, be dominant. In the case of F — bb, minimum-bias and other backgrounds 
may enter |7| . To the extent that the functional form (in the kinematical variables 
(p) of the sum of such backgrounds is substantially different from the /#(</>) cal- 
culated for the irreducible backgrounds, the sensitivities that we compute in this 
paper as a function of S/ B and S/ \f~B will have to be re-evaluated by the experi- 
mental collaborations. We re-emphasize that in the actual experiment /b(</>) will 
be directly measured using data for which the bb mass is outside the Higgs mass 
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peak region. 

In order to assess our ability to determine that a CP-even SM-like Higgs boson 
is indeed CP-even, we will follow the procedure of Ref. ]]J and consider three 
distinct Higgs coupling cases. We compare [recalling that a = a 2 /(a 2 + rb 2 )] 

• I) A Standard- Model-like Higgs boson, with a^O and 6 = 0, i.e. aj = 1, to: 

• II) a CP-mixed Higgs boson, with a = b, i.e. an = 1/(1 + r); and, 

• III) a pure CP-odd Higgs boson, with a = 0, b ^ 0, i.e. aui = 0, 

In order to compare with the results of Ref. HJ, we will compute the discrimination 
powers D12 = \olii,iii — aj\(S/yB)H(ai, Sj/B) defined in Eq. (H). These are a 
measure of our ability to ascertain that a CP-even SM-like Higgs boson is not CP- 
mixed or CP-odd, respectively. In comparing, we use only the best-case weighting 
functions Ocp considered in Ref. For F = 77 the best choice for Ocp was 
denoted by Oqp — &i; for F — bb the best choice was denoted by Oqp = b^. 

Table 1: ££77 channel, = 100 GeV: Discrimination powers D\ and D 2 for 
distinguishing a SM-like Higgs from a CP-mixed Higgs and from a CP-odd Higgs, 
respectively. Results from Ref. |1J are compared to those obtained using the optimal 
analysis. 





Technique 
b\ []J optimal 




Technique 
61 [0] optimal 


D 1 {B = 0)/y/S 

d 2 (b = oyVs 


0.137 0.376 
0.663 1.68 


D X (S = 130,5 = 21) 
D 2 {S = 130,5 = 21) 


1.47 3.68 
7.10 16.5 



Consider, first, F = 77 in the limit of B = 0. Table [I] gives D^ 2 (B = 0)/\^S, 
a first measure of the potential power for discrimination between the models; we 
see a large improvement relative to the Ocp = 61 weighting procedure by using 
the optimal techniques. Table || also shows the actual levels of discrimination that 
can be achieved if we adopt the LHC SM Higgs boson estimates of S = 130 and 
B = 21. Of course, since the total number of events available is not large, non- 
Gaussian effects in the statistics that are not incorporated in our approach will 
prevent achieving results quite this good. 

We turn now to the F = bb channel. The discrimination power D for separating 
models is determined as a function of S/ B using Eq. @ and H(ai, S/ B) (computed 
for the bb channel) as plotted in Fig. [l|. In Table |2], we give the number of signal 
events required to achieve D 12 = 2 for B/S = 50, and the corresponding S/\fB 
values. Also given are the results obtained using the best Ocp = 64 observable 
of Ref. []]]. The improvement using the optimal techniques is very remarkable; 
D 12 = 2 discrimination is possible with far fewer events and much lower S/ \f~B as 
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Table 2: tibb channel, nih = 100 GeV, B/S = 50: signal event rates, and cor- 
responding S/yB values, required to achieve Di j2 = 2 using the best Oqp = &4 
observable of Ref. |l| as compared to the optimal approach. 





D^B/S = 50) = 2 
S S/^fB 


£> 2 (£/S = 50) = 2 
S S/y/B 


h i 

optimal 


15000 17 
1592 5.64 


970 4.5 
80 1.26 




Figure 1: We plot H(a, S/B) for the ttbb final state as a function of S/B for the 
three cases: a = ai (solid); a = an (dashes); a = am (dots). 

compared to Ref. [T]. For large B/S, the results in the bb final state for D 12 are 
easily summarized: 

* = ", = 0,-*-, (10) 

with C\ = \au — ai\H = 0.36 and C*2 = \ani — ai\H = 1.61, where H = 1.61 is 
the result from Eq. (|9|). These values of C\ and C*2 can be compared to the values 
0.091 and 0.455 achieved in Ref. using Ocp = b^ the latter obviously imply 
much weaker discrimination power. 

A special remark is in order regarding the small S/\fB = 1.26 value required 
(see Table 0) for D 2 = 2 in the tibb channel using the optimal approach. If the 
Higgs boson mass is already known from observation in another channel (e.g. ^77), 
then one can perform the CP analysis for tibb events falling into the appropriate bb 
mass bin (assuming mass calibrations are well understood) even though a distinct 
bump is not observed. In this way, one could confirm and possibly improve CP 
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results obtained in the channel for which a mass bump is manifest. 

Although the 77 channel is not likely to be useful for Higgs bosons that do 
not have a large CP-even component, the bb channel is likely to be of very general 
utility. Light Higgs bosons will almost always have a large branching ratio for 
decays to bb. Given the large ttbb background rate, B/S will generally be large. 
Thus, it is useful to give the H(a, S/B) values, from which the discrimination power 
D of Eq. (El) can be computed, for other models. The results for H(au, S/B) and 
H(ciiii, S/B) appear along with those for H(ai, S/B) in Fig. |l|. At S/B = 1/50, 
all are of order > 1.55. They approach the common limit of H = 1.61. Using a 
lower bound of H(a, S/B) > 1.5 for B/S > 30 (the most likely range), we conclude 
that a marginal signal at the level of S/ \f~B = 3 will result in discrimination power 
of D > 4.5|a' — a\, regardless of the value a for the input model. This is quite an 
encouraging result in that for any Higgs boson that is observable in the ttbb channel 
one can hope to determine the mix of CP-even and CP-odd Higgs couplings to ti 
at a very useful level. 

Table 3: Ratio r = Ib/I a of cross section contributions arising from the CP-odd b 2 
term and CP-even a 2 term for various \y\ cut and rrih (in GeV) choices. 



Case 


\y\ < 4 

m h = 100 


\y\ <2 

m h = 100 


\y\ < 4 

rrih = 130 


M<2 
m h = 130 


r 


0.288 


0.243 


0.447 


0.362 



Table 4: ££77 channel: Di and D 2 discrimination powers obtained using the optimal 
technique, for the \y\ cut and rrih values considered in Table |[ 



Case 


\y\ < 4 


M<2 


\y\ < 4 


\y\ <2 




m h = 100 


m h = 100 


rrih = 130 


rrih = 130 


D^B = 0)/v^S 


0.376 


0.292 


0.449 


0.336 


D^S = 130,5 = 21) 


3.68 


3.07 


4.31 


3.50 


D 2 (B = 0)/y/S 


1.68 


1.49 


1.45 


1.26 


D 2 (S = 130,5 = 21) 


16.5 


15.7 


14.0 


13.2 



We now explore sensitivity of these results to the \y\ cut and to m h . Since 
the \y\ < 4 rapidity cut for the t, t and h may be too generous to allow for full 
reconstruction in the relatively complicated tth final state, we also consider \y\ < 2. 
We will explore sensitivity to rrih by giving results for rrih — 130 GeV. First, we 
give the r = l h j ' I a values, appropriate in all the cases considered, in Table |^. Note 
that as rrih increases the overall signal cross section contribution (/{,) of the CP-odd 
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b 2 term increases relative to the contribution (J ) of the CP-even a 2 term, implying 
larger r values at larger m^. Next, in Table | we give the discrimination powers 
D12 in the 77 mode. For fixed values of S = 130 and B = 21, the stronger \y\ < 2 
cut gives D12 values that are not so much smaller than for \y\ < 4. However, this 
masks the fact that only 1/3 as many events are accepted for \y\ < 2 as compared 
to \y\ < 4. Reducing S and B by a factor of 3 (for \y\ < 2) from S = 130 and 
B = 21 would yield D 12 values that are a factor of l/y/3 smaller than given in 
Table |j. Regarding nih dependence, note that, keeping S and B constant, D\ (D 2 ) 
is larger (smaller) for nih = 130 GeV as compared to nih = 100 GeV. The larger 
D\ is because \otu — oti\ is larger (larger r) at nih = 130 GeV. This contrasts with 
the case of D2 for which \aju — ai\ = 1 independent of ro/,. 

Table 5: ttbb channel: signal event rates, and corresponding S/ yB values, required 
to achieve Di 2 = 2 using the optimal approach. 





D 1 (B/S 


= 50) = 2 


D 2 (B/S 


= 50) = 2 


(\y\ Cut,m h ) 


S 


s/Vb 


S 


s/Vb 


(4,100) 


1592 


5.64 


80 


1.26 


(2,100) 


3215 


8.02 


123 


1.57 


(4,130) 


932 


4.32 


89 


1.33 


(2,130) 


2018 


6.35 


143 


1.69 



Table 6: ttbb channel: asymptotic H values and corresponding values of C% and 
C 2 , as defined in Eq. fllUp. 



Case 


\y\ < 4 

m h = 100 


\y\ < 2 

m h = 100 


\y\ < 4 

nih = 130 


M<2 
nih = 130 


H,C 2 
Ci 


1.61 
0.360 


1.29 
0.252 


1.53 
0.473 


1.20 
0.319 



Some comparative results in the ttbb channel are given in Tables [| and [!]. The 
H(a, S/B) values for nih = 100 GeV and \y\ < 2 were given in Fig. [p. For S/B = 
1/50, one finds H(aj, 1/50) = 1.27, H{a n , 1/50) = 1.27, and H(a IH , 1/50) = 1.24, 
a significant deterioration from the > 1.55 values found for \y\ < 4. The S/B — > 
limit is H = 1.29, again significantly smaller than the 1.61 value obtained with 
\y\ < 4 cuts. More limited acceptance implies that more signal rate (which is 
harder to achieve with more limited acceptance, e.g. \y\ < 2 accepts only 1/3 
as many events as \y\ < 4) and larger S/y/B is required for good discrimination 
between the various models. Clearly it is important to accept as much of the phase 
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space as possible, both to maximize the signal rate and to maximize the region over 
which shape differences between models can be explored. The trends with are 
similar to the 77 channel, with Di (D 2 ) being larger (smaller) for = 130 GeV 
as compared to rrih = 100 GeV, for the same S/ \f~B. 

3 The MSSM Higgs Bosons 

A Higgs boson sector of particular interest is that of the minimal supersymmetric 
model (MSSM) §. If the mass of the CP-odd A of the model is > 150 GeV, then 
the light CP-even h° will be SM-like, and our tth analysis will apply in both the 
77 and bb channels. However, the heavier CP-even H° and the A will probably 
present a greater challenge. If the tan /3 parameter of the MSSM is near 1, then the 
tiH ,tiA° couplings, and corresponding LHC production rates, will be substantial 
for modest m H o,m^o. Depending upon how easily they can be discovered in the 
relevant decay channels (such as H°, A — > bb, t + t~ , H° — * h°h°, A — > h°Z, and, 
for m#o, m^o > 2m t , H°, A — > tt) and upon how small S/B in each channel is, 
our tth procedures could be useful. We do not attempt a detailed study here. 

If the tan/5 parameter of the model is large and m^o ~ m H o > 150 GeV, then 
the tiH and ttA° couplings are suppressed by a factor of 1/ tan (3 and those to 
bb (and r + r~) are enhanced by a factor of tan/3. As a result, the primary LHC 
production modes are pp — > bbH°X, bbA°X, where both H° and A decay to bb 
(t + t~) with a branching fraction of order 0.9 (0.1). In this case, only the bbbb 
and bbr + T~ final states are relevant for detecting the H° and A . Both ATLAS 
H and CMS || claim that S/\/B = 5 discovery of the H°,A° is possible (for full 
'three-year' L = 300 fb _1 integrated luminosity per detector) in the bbr + r~ mode 
provided tan/3 is above a m^o-dependent lower limit ranging from tan/? ~ 2 at 
m^o ~ rriHO ~ 150 GeV up to tan (3 ~ 20 at m^o ~ m#o ~ 500 GeV. Prospects 
in the bbbb final state are more controversial. The bbr + r~ mode is preferable for 
our purposes in any case because the r's allow us to identify the Higgs boson 
without any combinatoric uncertainty. The question, then, is whether or not our 
techniques will allow verification of the CP= + and CP= — nature of the H° and 
A , respectively. 

First, and very crucially, we do not need to identify b vs. b due to the fact 
(see earlier) that f a ((fi) and fb((fi) are independent under interchange of the b and b 
momenta. However, there are two significant difficulties. First, in Eq. (|I|) / a (0) ~ 
fb{4>) in the limit where the quark mass approaches zero. (That is, S(0) oc a 2 + b 2 
in this limit, the a 2 — b 2 term being suppressed by ~ "^quark/v^i in comparison.) 
Consequently, a very large number of events will be required for a statistically 
useful determination of the a 2 — b 2 coefficient. Second, mass resolution in the t + t~ 
(or bb) channel will typically not be adequate for separating the H° and A mass 
peaks, implying that one could at best verify that a 2 — b 2 = after averaging over 

3 Their 77 branching ratios are typically too small to be useful. 
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the H° and A ( given that the H° and A have the same production rate at large 
tan/3). Our first estimates suggest that, for tan/5 values such that discovery is 
possible, a statistically meaningful verification of the a 2 = b 2 expectation relative 
to pure a 2 = or pure b 2 = might be possible, despite the small m&/ vl coefficient 
multiplying a 2 — b 2 . A careful study of this issue by the detector collaborations - 
i.e. one that includes all resolutions, minimum-bias backgrounds, and so forth - 
is needed to properly assess the prospects. 

4 Final Remarks and Conclusion 

In this paper, we have employed the optimal analysis procedures of Ref. |2j to 
determine the accuracy with which the relative magnitude of the CP-even and 
CP-odd it couplings of a light Higgs boson can be determined using the final 
state kinematical distributions of the t, i and h relative to one another in LHC 
pp — > tth events. We have found that the optimal analysis procedure yields a large 
improvement in the statistical discrimination power (D = v5?) as compared to 
the simpler weighting function techniques previously considered in Ref. fjj. 

We emphasize that we have employed a discrimination measure that does not 
make use of the (substantial) sensitivity of the absolute production rate to the 
CP-even versus CP-odd couplings; this is desirable since absolute rates have ad- 
ditional systematic uncertainties associated with model dependence for branching 
ratios, experimental efficiencies and the like. However, an estimate or, preferably, 
experimental measurement of the signal to background ratio is needed in order to 
compute quantitative results for the x 2 associated with a model that is different 
from the input or true model. 

It is useful to note that in the actual experimental analysis an equivalent of the 
optimal analysis procedure employed here would be to determine the likelihood 
of the fit, as a function of the a and b parameters, to the detailed phase space 
distribution (allowing the overall normalization to be a free parameter). In the limit 
of Gaussian statistics, this and the theoretically much more convenient optimal 
procedure employed here are equivalent: the ratio of the likelihood of a fit in an 
incorrect model to that in the correct model will be equal to exp[— x 2 ], where the 
X 2 is the x 2 °f an incorrect model choice as computed in the optimal procedure 
using inputs based on the correct model. 

Our conclusion is that if a Higgs signal is observable in the tth (h — > F, F = bb 
or 77) final states at the LHC, then the relative size of the CP-even and CP-odd 
couplings of the h to it can be determined at a statistically useful level. In particu- 
lar, in the case of a SM-like Higgs boson, we find that Higgs CP coupling mixtures 
that are significantly different from pure CP-even can generally be excluded at a 
high level of statistical significance. 

Finally, we noted that there is reason to hope that the techniques studied here 
could be employed for studying the CP nature of the heavier Higgs bosons of 
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the minimal supersymmetric model. In particular, for tan (3 values large enough 
that a signal is observable in the bbr + r~ final state coming from bbH° and bb~A° 
production (with H°,A° — > t + t~), one can hope to verify that the overlapping 
H° and A resonances yield the expected equal mixture of CP-even and CP-odd 
signals. 
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